Self-consistent model of ultracold atomic collisions and Feshbach resonances in tight 
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We consider tlie problem of cold atomic collisions in tight traps, where the absolute scattering 
length may be larger than the trap size. As long as the size of the trap ground state is larger 
than a characteristic length of the van der Waals potential, the energy eigenvalues can be computed 
self-consistently from the scattering amplitude for untrapped atoms. By comparing with the exact 
numerical eigenvalues of the trapping plus interatomic potentials, we verify that our model gives 
accurate eigenvalues up to milliKelvin energies for single channel s-wave scattering of ^"^Na atoms in 
an isotropic harmonic trap, even when outside the Wigner threshold regime. Our model works also 
for multi-channel scattering, where the scattering length can be made large due to a magnetically 
tunable Feshbach resonance. 
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I. INTRODUCTION 

Along with the development of laser cooling of atoms 
have come techniques for trapping the cold atoms, with 
tremendous advantages for experimental atomic physics. 
Just to name a few potent examples where trapping is 
necessary, the Bose-Einstein transition has been reached 
in several atomic species [1 J3|, threshold scattering prop- 
erties have been studied [ i pT, molecules formed with the 
assistance of light ^ , and quantum chaos and quan- 
tum phase transitions observed using optical lattices, 
a 

Atomic collisions play an essential role in most of these 
phenomena. In the past one could ignore the fact that 
these collisions take place in a trap, since trap sizes 
are very large in comparison with the sizes associated 
with atomic interactions. However, recent developments 
make it crucial to account for the effect of trap confine- 
ment on collisions when the atoms are held tightly in 
one, two, or three dimensions by optical lattices. For 
example, Greiner et al. Q have observed a quantum 
phase transition from a superfluid to a Mott insula- 
tor within a three-dimensional optical lattice. More- 
over, several low-dimensional phase transitions of cold 
bosonic systems have been conjectured. In two dimen- 
sions, the Kosterlitz-Thouless transition may occur [ p^ , 
while in one dimension the Tonks-Girardeau phase should 
be possible [|ll|jl2[ . Zero-temperature transitions have 
also been investigated theoretically for a rotating two di- 
mensional gas |13| . All of these transitions depend on 
atomic collisions, and for quantitative predictions the 
low-dimensional interactions must be understood. 

Two proposals for quantum computing involve loading 
cold atoms into optical lattices, and using the interaction 
between the atoms as the switching mechanism p^ , p^ . 
In one type of quantum logic gate, two atoms are brought 
together and allowed to interact for a set time interval, 
resulting in different phase shifts depending on their hy- 



perfine sublevels. A recent experiment Q represents an 
important first step towards quantum logic applications, 
since it shows that a lattice can be initialized with uni- 
form occupancy of lattice sites. 

Another burgeoning area is the study of Feshbach res- 
onances, and weakly-bound molecular states, in the in- 
teraction of two ultra-cold atoms. This has improved the 
knowledge of interaction parameters of alkali atoms and 
opened up the field of molecular condensates and 3-body 
processes jl^,^ . By tuning Feshbach resonances one can 
easily reach an interesting regime where the scale length 
associated with the scattered wave exceeds the trap width 

We address these problems by calculating the eigen- 
values of two interacting atoms confined in a trapping 
potential. A popular method for representing cold atom 
interactions is to replace the exact interatomic potential 
by a delta-function pscudopotential proportional to the 
scattering length a: 



V = a5[Y)—r, 

m or 



(1) 



where m is the atomic mass and r is the interatomic sep- 
aration. An analytic solution for the eigenvalues of two 
atoms in an isotropic harmonic trap plus the pscudopo- 
tential Eq. (|l|) has been found However, some of us 
have previously shown that the use of this solution is lim- 
ited to sufficiently weak traps such that the trap width is 
much larger than \a\ Here we reexamine this prob- 
lem and propose a self-consistent method of calculating 
the trap energies, which gives good quantitative results 
over a wide range of trap frequencies, even when \a\ is 
larger than the trap size. The essence of our model is 
to replace a with an energy-dependent effective scatter- 
ing length. An advantage of our model is that once the 
energy-dependent scattering phase shift for a particular 
type of cold collision is known, either from experiment 
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or from close-coupling calculations, it can be easily ap- 
plied to obtain eigenvalues for traps of all frequencies. 
Conversely, if the eigenvalues are measured, information 
about collisions can be obtained. 

We note that the pseudopotential can be used to ob- 
tain approximate solutions for trapped colliding atoms in 
one dimension [ pT|JTl| and two dimensions p^ j. It may 
be possible to adapt our self-consistent method to accu- 
rately treat scattering in "cigar-" or "pancake-" shaped 
traps. 

The paper is organized as follows. In Sec. II we for- 
mulate the problem of atoms colliding in a tight spher- 
ical trap, and briefly review scattering theory. In Sec. 
Ill we motivate and explain the self-consistent eigen- 
value model, which is our main result. Limitations of 
the model are discussed. Section IV applies the model 
to single-channel scattering of ^^Na atoms in a trap, and 
shows good agreement with numerical calculations using 
the full interaction Hamiltonian. Section V demonstrates 
similar good agreement for the case of multi-channel scat- 
tering. Specifically, we consider a magnetically-induced 
Feshbach resonance in Na2. We compare, for a range 
of magnetic fields, exact numerical results from the five- 
channel close-coupled scattering problem with the self- 
consistent eigenvalues. Finally, in Sec. VI we draw con- 
clusions and consider more general traps and applications 
to many-body theory. 



II. TWO ATOMS COLLIDING IN AN 
ISOTROPIC HARMONIC TRAP 

We consider an isotropic harmonic trap described for 
atom J = 1, 2 at position Vj by 



Hrapl^j) 



1 



2 2 

mcu r- , 



(2) 



where lo is the trapping frequency. Harmonic traps can 
be made by a variety of means. Very tight confinement 
is possible with a three-dimensional optical lattice. Typ- 
ical experimental trap frequencies range from 50 kHz to 
1 MHz. These optical dipole traps are much tighter than 
those obtained with magnetic fields. In a recent experi- 
ment 1^, isotropic potentials at each site were produced 
from three optical standing waves of equal intensity. 

For the isotropic harmonic trap, the two-atom Hamil- 
tonian is separable in the center-of-mass and relative co- 
ordinates. Since the center-of-mass motion is just that 
of the well-known isotropic harmonic oscillator, we need 
only discuss the problem in the relative coordinates. The 
Hamiltonian is 



H = 

2^ 



^M^'^' + ^intW, 



(3) 



where r = |ri — r2 1, = m/2 is the reduced mass, and 
Vjj^^(r) is the interatomic potential. In (relative) spher- 
ical coordinates, the trap states neglecting V\^if^{r) have 
energy eigenvalues 



E'^?^ = ( 271 + L + ^ ) 



(4) 



where n = 0, 1, 2, ... is the radial quantum number and 
L = 0, 1, 2, ... is the partial wave quantum number. We 
henceforth consider only s- waves [L = 0). The size of 
the ground state trap wavefunction is characterized by 



(5) 



Typical trap sizes I for Na in the above mentioned trap 
frequency range are 30 nm to 130 nm. 

The interatomic potential Vy^i^{r) is characterized by 
a short-range region of strong chemical bonding and a 
long-range van der Waals potential. 
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and leads to a van der Waals scale length, P,p3U2 



1 /2/iC6 



"° = 2 



1/4 



(6) 



(7) 



For r <^ xq the scattering wavefunction oscillates rapidly 
due to the strong interaction potential. In alkali ground 
state interactions, Cg is the same for all hyperfine states 
of a given atomic pair; consequently, is the same for all 
collision channels. In the case of Na2 considered below, 
it is about 2.4 nm. 

For collisions of atoms in the absence of a trapping po- 
tential, the asymptotic s-wave scattering wavefunction 
for relative collision momentum fik approaches 



sin(fcr -f 6q) 
\fkr 



(8) 



at large interatomic separation r xq. Another length 
scale that naturally appears for cold collisions is the scat- 
tering length, defined in terms of the s-wave phase shift 
<^o by 



lim 

fc-*0 



tan5o(fc) 



(9) 



The Wigner law regime is then defined by the range of 
momenta for which = ka is a good approximation, i. 



e., 



k < 



2\a\ 



(10) 



The scattering length can take on any value between +oo 
and —oo. As \a\ becomes large, the range of k for which 
the Wigner law applies becomes very small. 

In view of typical trap sizes and van der Waals length 
scales, we need only consider the experimentally accessi- 
ble regime, for which 
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xq <C L 



(11) 



On the other hand, the scattering length can have a larger 
magnitude than the trap width I. This is especially likely 
if the scattering length is modified by means of a Fesh- 
bach resonance. 

Our goal is to find a simple model for calculating the 
new energy eigenvalues of the trap when collisions are 
present. An analytical solution of this problem was pre- 
sented in Ref. l2C|l by replacing V^^.^ by the pseudopo- 



tential of Eq, 
Wigner law is va. 



This replacement assumes that the 
lid. However, we previously showed that 
the eigenvalues thus obtained are not always in agree- 
ment with numerical results | [ig[ |! Specifically, they are 
least accurate when \a\ approaches or exceeds I. One way 
to see this is that the energies of the unperturbed trap 
states are already large enough that the Wigner thresh- 
old law is invalid. For the unperturbed trap ground state 
E = 3huj/2 and hence the root- mean-square momentum 
k = VS/l. Therefore by Eq. M) we are outside of the 
Wigner regime if \a\ > 7r/(2V3)L 

In the next section we will use the inequality Eq. 
to motivate a self-consistent model of cold collisions in 
the trap, that is valid at all relevant energies and scat- 
tering lengths. 



III. SELF-CONSISTENT MODEL 

The improved model we propose relies on a general- 
ization of the pseudopotential approximation for V^.^^ in 
Eq. (p. We introduce the energy-dependent pseudopo- 
tential operator 129] 



^eff = ^«eff(^)^W^^' 
where the effective scattering length is defined as 



«eff(^) 



tan(5o(A:) 
k ■ 



(12) 



(13) 



and the kinetic energy is related to the momentum by 
E = IV^k^ /2fi. This operator gives the same asymptotic 
wavefunction, Eq. (^, as the full interaction potential 
V^jj^ . The effective scattering length reduces to the usual 
one, Eq. (H), in the Wigner threshold regime. The phase 
shift in Eq. ( p^ does not need to be small in order to 
use Eq. ([l^) . Even though the effective scattering length 
diverges when Sq is an odd multiple of 7r/2, the wave- 
function remains well-behaved. 

Reference found the eigenvalues of the trapped 
atoms interacting through the operator in Eq. (jl|) as the 
solutions of the equation, 



y = /(n 

where the "intercept function" is 



(14) 



, , 1 fnE ^\ r(J- + i) 



(15) 



and r is the gamma function. To account properly for the 
scattering in tight traps, where the Wigner law may not 
apply at the trap energies, we need to replace Eq. (|lj) 
by one in which the left-hand side is energy-dependent 
and solve the equation 



«eff(^) 



I 



f{E) 



(16) 



self-consistently for the eigenvalues. 

One might ask, why does the idea of the pseudopo- 
tential still work outside the regime of the Wigner law? 
The answer is that the collision occurs on the very short 
length scale ccq, so the interatomic interaction potential 
is undistorted by the trap. This in turn means that the 
kinetic energy at which the effective scattering length 
needs to be evaluated is the eigenvalue itself, since the 
trap potential is negligible for r < xq. Thus we were led 
to the self-consistent Eq. (p^). 

This model can be expected to break down if the trap 
becomes too tight. The interatomic potential V-y^i^ be- 
comes comparable to the trap potential near r = ^/Txq. 
Hence the inner part of the wavefunction where the scat- 
tering occurs is nearly the same as that without the trap 
when Xq ^ \/lxo, equivalent to Eq. ( pi] ) . A different kind 
of limitation is that this model cannot predict bound 
states without our knowing the analytical continuation 
of the effective scattering length to negative energies. 



IV. SINGLE-CHANNEL SCATTERING 

The first problem we consider is that of doubly polar- 
ized (electron and nuclear spin up) ^^Na atoms colliding 
in the trap. In this case, there is only one scattering chan- 
nel, governed by the a'^S+ adiabatic Born-Oppenheimer 
potential. The scattering length is a = 3.2 nm, and Figs. 
|l| and ^ show the effective scattering length as a function 
of energy. It increases with energy and diverges near 
E/h = 90 MHz where Sq = 7r/2 (this corresponds to a 
local maximum of the s-wave cross section), and is neg- 
ative immediately above this energy. In this work both 
the single- and multi-channel phase shifts are calculated 
by applying the Gordon propagation method [E6| with 
the best available scattering potentials for Na2 [P7| . 

The radial Schrodinger equation for the Hamiltonian 
Eq. (^ was solved numerically for the eigenvalues. For 
a detailed description of our numerical method using a 
discrete variable representation, see Ref. |2^. We take 
a trap frequency of uj/2n = 1 MHz, for which I — 29.6 
nm and Tito /ks = 48 /LtK [ks is the Boltzmann constant). 
Such a tight trap should be feasible in a Na optical lat- 
tice. 

We illustrate the graphical solution of the self- 
consistent model in Figs. |l| and ^. In each plot, the solid 
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curve represents the left hand side of Eq. (^^, 
while the dashed curve is the right hand side. The ab- 
scissae of the points where the curves intersect give the 
self-consistent eigenvalues according to the model. One 
way of comparing with the exact numerical eigenvalues 
is to evaluate the intercept function f{E) at these eigen- 
values; these points are plotted as circles. The closer 
the circles lie to the intersection of the curves, the better 
the agreement. The exact numerical and model eigenval- 
ues in Fig. |l| agree to better than 0.0016 Tiw. Note that 
the solution of Eq. (|lj) is found from the intersection of 
f{E) and the horizontal line E = a/l; the corresponding 
eigenvalues differ significantly from both the exact and 
self-consistent ones. 
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FIG. 1. Effective scattering length (solid curve) and inter- 
cept function f{E) (dashed) versus energy for doubly polar- 
ized ^^Na in a 1 MHz trap. The energies at which the two 
curves intersect give the self-consistent eigenvalues. The cir- 
cles show the actual positions of the exact numerical eigen- 
values along the curve of the intercept function. 



The range of energies in Fig. is centered near the en- 
ergy at which a^g/l diverges. Even though la^g] I, the 
self-consistent eigenvalues are still accurate. They agree 
with the exact ones to within 0.0018 Tiuj. Clearly this val- 
idates our model. We have also obtained eigenvalues for 
much higher trap frequencies, at which distortion of the 
collision potential is expected to cause the self-consistent 
model to fail. At a trapping frequency 100 MHz, where 
I = 2.96 nm « xq , the error between the exact eigen- 
values and those obtained from our model has increased 
to 0.045 hcu. The crucial interaction length scale for com- 
parison to the trap size I is xq, not the effective scattering 
length Ogg. 
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FIG. 2. Same as Fig. |l| but at a higher energy range. 
Note that the effective scattering length diverges, but the 
self-consistent eigenvalues still agree with the exact numer- 
ical eigenvalues. 

The difference between the lowest seven eigenvalues 
and the corresponding harmonic oscillator eigenvalues 
given by Eq. (^ is plotted in Fig. ^ versus the quantum 
number n. The shift due to the interactions is a signif- 
icant fraction of huj and should be observable in appro- 
priate experiments. The dependence of the shift on the 
index for the lowest few eigenvalues is due mainly to the 
energy dependence of the gamma functions in Eq. (|T^), 
and only slightly due to the variation of the effective scat- 
tering length with energy. On the other hand, for the 
higher eigenvalues in Fig. ||, the shifts in eigenvalues 
arise mostly from the rapid variation of effective scatter- 
ing length with energy. Near the asymptote a^Q oo 
the eigenvalues have increased by approximately huj com- 
pared with the unperturbed values. 

The above examples show that accurate eigenvalues 
can be obtained by using results of the single-channel 
scattering problem (without the trap), and solving the 
self-consistent Eq. (iq). Our self-consistent model is good 
even when the effective scattering length is large com- 
pared to the trap width, provided the trap size is still 
larger than the van der Waals length scale. 



MULTI-CHANNEL SCATTERING AND 
FESHBACH RESONANCE 



In the previous Section, large ratios of effective scatter- 
ing length to trap width were only possible for very high- 
lying levels. Here we want to discuss a situation where 
jflggl/^ is arbitrarily large for the lowest trap levels. This 
can be experimentally realized for s-wave collisions using 
a magnetically-tuned Feshbach resonance. 
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effective scattering lengtfi Eq. (|T^ for the {aa} cliannel 

is 



2 3 4 

Quantum number n 

FIG. 3. Difference between eigenvalues for interacting and 
noninteracting, doubly polarized Na atoms in a 1 MHz trap 
versus quantum number n. The self-consistent and exact nu- 
merical eigenvalues are indistinguishable on the scale of this 
figure. 

We consider a Feshbach resonance in the collision of 
two ^^Na atoms in their lowest hyperfine level at a mag- 



netic field near 90 mT [|17|,^ 30 1 . The hyperfine states of 
the ^'^Na atom diagonahze the Zeeman and hyperfine in- 
teraction and are labeled by \a) , \b) , . . . \ h) , starting from 
the lowest internal energy. For very low collision en- 
ergy, s-wave collisions of two \a) atoms are represented 
by five symmetrized asymptotic collision channels, one 
of which is open, |{aa}), and four of which are closed, 
|{a5}), \{fh}), and \{gg})- The interaction be- 

tween the atoms is mediated by the X^I]+ and a'^S+ adia- 
batic Born-Oppenheimer potentials. During the collision 
this interaction mixes hyperfine states and is described 
by a Hamiltonian coupling the above five channels 
A Feshbach resonance state at energy Ep is located at 
the threshold of the |{aa}) channel for a magnetic field 
Bxes ~ 90.09 mT. This resonance is a quasibound molec- 
ular eigenstate of the four closed channel problem. It can 
be formed from or decay to the |{aa}) open channel, to 
which it is coupled. As the magnetic field B is changed 
near i?res: the resonance energy also varies with B: 



The analytic theory of Feshbach resonances [^JSl) 
shows that the phase shift Sq can be written as the sum 
of background and resonant scattering contributions: 



-^0 



(5|-,g — arctan 



r, 



2{E - Ep - Ap)' 



(18) 



where F^? is the linewidth, A^? is a level shift induced by 
the coupling between the open and closed channels, and 
(5]-,g is the background phase shift. It follows that the 



«eff(^) 



^ - {E - Ep - Ap)tanS^g 
k {e - Ep - Ap + '^i&n 5^^ 



(19) 



Up to the highest energy we will consider, E/h f^i 5 MHz, 
both Fp and tan(5|-,g are proportional to -s/S, and Ap 

becomes constant. Moreover, Eq. (|l^) shows that the 
effective scattering length diverges near the energy 



-^div - 



Ep + Ap 



1 



1 d 

2 dE 



(r. 



tan 6] 



bg 



(20) 



E^O 



The effective scattering length is positive below and nega- 
tive above -E^jj^, which is magnetically tunable according 
to Eq. (17). However, instead of employing the analytic 



theory, at a given value of magnetic field we directly ob- 
tain the effective scattering length as a function of E from 
a numerical close coupled scattering calculation with five 
channels. This enables us to extract the position of the 
divergence E^^^, which is plotted as the dashed curve in 
the {E, B) plane in Fig. |. 
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FIG. 4. Numerical (circles) and self-consistent (solid 
curve) eigenvalues vs magnetic field B for Na in a 500 kHz 
trap. The dashed line shows iS^jj^, where the effective scat- 



(17) tering length diverges at a fixed value of magnetic field. 



We now examine the effect of the Feshbach resonance 
on trap eigenstates, assuming a trap with lo/2'k = 500 
kHz. We used the numerical discrete variable method 
for five channels to calculate the lowest eigenvalues of the 
trap states for a range of magnetic fields near the reso- 
nance. These eigenvalues are plotted as circles in Fig. IJ. 
Solutions to the self-consistent eigenvalues were obtained 
by solving Eq. (16) graphically as in the single-channel 
case; the solutions versus magnetic field are the solid 
curves. The self-consistent eigenvalues agree well with 



5 



the numerical ones for all values of energy and magnetic 
field; the worst agreement, < O.lhuj, is for eigenvalues 
near the resonance position. The self-consistent eigen- 
values always lie slightly above the numerical ones. Note 
that the self-consistent eigenvalues cross the E^^^ curve 
near E/Tilo — 1/2, 5/2, 9/2 . . . Another particular feature 
of the plot, which is correctly reproduced by the numer- 
ical solution, is that as B decreases the lowest trap state 
[E > 0) becomes the highest bound state {E < 0) for 
a magnetic field B < iJres- This occurs where the ef- 
fective scattering length is still finite and positive, since 
a^Q/l « 1.48 when E = Qm Eqs. 
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VI. CONCLUSION 

We have shown how a self-consistent model can be used 
to calculate the eigenvalues of interacting atoms in an 
isotropic harmonic trap. Our model involves solving an 
equation containing the effective scattering length for un- 
trapped atoms, and the trap frequency. We compared our 
model with exact results for ^^Na both for a single chan- 
nel collision and a multi-channel collision with a tunable 
Feshbach resonance. In both cases, the model can accu- 
rately treat tight traps, as long as the trap size is larger 
than the van der Waals scale length. Consequently we 
expect the model to apply to other atomic species. In 
particular, Cs would be an interesting case for which the 
scattering length is large in comparison with even modest 
trap sizes |3^ . 

In the future, we want to generalize the self-consistent 
model to more arbitrary trap potentials. There are two 
technical problems to be overcome. First, for atomic col- 
lisions in anisotropic harmonic traps, the relative coordi- 
nate equation does not separate; this implies that differ- 
ent partial waves are coupled via the anisotropy. A re- 
lated point is that the scattering of higher partial waves 
can also be modeled by pseudopotentials Second, for 
anharmonic traps, the center-of-mass and relative atomic 
coordinates do not separate, and even more coordinates 
must be treated simultaneously. Anharmonic terms be- 
come important for low intensity optical lattices or for 
trap levels with many quanta of excitation. 

One would expect to be able to use the effective scat- 
tering length in many-body problems, where the pseu- 
dopotential approximation has had widespread use. One 
would simply need to replace a by the effective scattering 
length. This should be especially useful and necessary for 
situations where a tunable Feshbach resonance is used to 
alter the interaction properties. There are a number of 
cases where the relative collision energy for a many-body 
system is well-defined, such as for condensates in optical 
lattices colHding condensates p3| , ^ , or cold gases 
of mixed fermionic species, where collisions occur at the 
Fermi energy. It should also be possible to incorporate 
inelastic collision loss channels by using a complex effec- 
tive scattering length [Q. 
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